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properties of spectral fluctuations of quantum-chaotic systems is considered 
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I. INTRODUCTION 



A quantum graph system consists of a quantum particle moving along the bonds of an 
arbitrary finite graph G [1]. In the classical limit, this system generates a simple stochastic 
dynamics, which is specified by the translational motion along the bonds of the graph and 
stochastic scattering at its vertices with preset scattering probabilities. This dynamics has 
many common features with the dynamics of usual chaotic systems [2]. For example, periòdic 
trajectòries in such a system are isolated and their number increases exponentially with the 
period. At the same time, the statistical behavior of various spectral characteristics of 
sufficiently complex quantum graphs, e.g. the probability distribution of spacings s n = 
k n — k n -i between the nearest levels of the momentum was numerically shown [1] to follow 
the predictions of the Random Matrix Theory (RMT) [3,4], as it is usually the case for 
classically nonintegrable systems. It also turns out that a great number of problems of 
classical and quantum dynamics on the graph allow exact solutions, which makes these 
systems convenient models in the context of the analytical theory of "quantum chaos". In 
particular, for these systems there exist the exact periòdic orbit expansions of the quantum 
density of states (Gutzwiller formula) [1] along with a similar expansion for the spectral 
staircase: 

N(k) = J2 Q ( k ~ k j) = N(k) + - Im A p e iL ? )k , (1) 
j=i " p 

Here N(k) is the average number of levels in the range [0, k], L p ^ is the optical length of the 
periòdic trajectory with the index p, and A p is a certain weight factor explicitly defined in 
terms of the scattering coemcients at the graph vertices. It should be emphasized that the 
existence of the explicit expansions of the global characteristics such as (1) is not equivalent 
to the ultimate solution of the spectral problem, which should provide local information 
about the individual levels in the form of an explicit dependence k n = k{n). An approach 
for determining the quantities k n explicitly was proposed in [5], which is based on using a 
finite system of r + 2 auxiliary "separators" kn \ kn+v ■■■■> k n +1 \ the first of which is the 
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physical spectral sequence k n — k n , and the last one is a globally defined explícit function 
of n: 

The key property of these sequences is that they must satisfy the "bootstrapping" conditions 

%{í) <- <- i.0') (qN 

which guarantee that between every pair of the neighboring points and fcjf2i (see Fig. 1) 
there exists a single point kn~ l \ In [5] it was also pointed out that due to certain analytical 
properties of the spectral determinant A(k) = 1 + 'Yl li aie tkL( - i \ where Lu), are different 
linear combinations for the bond lengths í a , Z 2 , In b i the set kn can be provided by the 
sequence of zeros of the j-th derivative of the function A(k) [5,8]. In this case, the quantity r 
characterizing the degree of spectral irregularity is defined as the minimal number for which 
the condition |aj (L(j)/L ) r | < 1 is satisfied [5]. 




FIG. 1. Bootstrapping of spectral staircases for separating sequences kn of the completely 
connected four- vèrtex graph with r = 7. The plots N^\k) are vertically shifted for the sake of 
clarity. It is clear that the physical spectral staircase N^°\k) is interlaced by the staircase N^'Qe), 
etc. The last staircase N^ r \k) is intersected by the Weyl average N(k) 

In the simplest case of regular graphs when r = [6,7], only one auxiliary sequence (2) 
is required and various spectral characteristics can be calculated using the formula 
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f(k n )= f(k)p(k)dk. (4) 

As pointed out in [6,5], this case corresponds to the situation in which the straight line 
with the slope Lo/n, representing the Weyl average N(k), "pierces" the physical spectral 
staircase N(k), i.e., N(k) intersects every stair step of N(k) at the points kn\ 



II. STATISTICAL PROPERTIES OF THE SPECTRA OF REGULAR GRAPHS 

Using the Gutzwiller formula in Eq. (4), one can derive the explicit expansions for 
various spectral characteristics fn\ for example, for fructuations 5^ = ^ (k n — k n ), of the 
eigenvalues k n around the Weyl average or for the distances between levels s„ im = k n+m — k n . 
Such expansions have the form [6,7] 

/ (o) = / (o)_^ C7 (o) cos( ^o) n + V9 (o) )) (5) 

V 

where the freqüències iJf 1 are defined via the periòdic orbit lengths as, uf^ = 
The first term of expansion (5) determines the average value of the quantity fn \ whereas 
the following sum describes fluctuations around the average. Each frequency oo^ is an 
integer combination uj^ = m^Jíli + mj^^ + ••• + m ^\í B ^N B , of the quantities ílj, which 
are expressed in terms of the lengths of the graph bonds as íl, = k/Lo, and the coefficients 
rrip - indicate how many times the orbit passes along the bond k- The sum 
m p°2 + ••• + m p°Af B -i specifies the total number of scattering events that the particle moving 
along the trajectory p undergoes at the vertices. If Eq. (5) includes only the orbits for which 
\mf\ < m, we arrive at the m-th approximation to the exact value fn^ [1,6]. 

Since the numbers íl, satisfy the condition Q 1 + íl 2 + ... + £In b = 1, only N B — 1 of these 
numbers are independent. Expressing one of them, e.g., fÍN B , in terms of the others, let us 
consider the (genèric) case when the numbers ílj = Qí — Q Nb are irrational and algebraically 
independent. Let us call the orbit p algebraically simple (with the notation p') if the integer 
coefficients mi°] = mi°] — mi°L have no common divisors. Such orbits in general differ from 
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the dynamically simple orbits that correspond to single traversals along closed sequences of 
bonds during the partides motion along the graph [1,2,5-7]. 

The expansion (5) enables one to pass immediately to the statistical description of the 
sequence fn \ Indeed, it is well known that the sequence of the remainders x n = [cm] modl 
for any irrational number a and n — 1,2, is uniformly distributed in the interval [0, 1] 
[14]. Since the arguments of the trigonometric functions appearing in series (5) are defined 
modulo 2n, parsing through the vàlues fn^ yields a sequence which is statistically equivalent 
to the series 

/i 0) = / (0) -E^ 0) ^(mfx + v,f), (6) 
v 

Here, Òf^ and rh^ correspond to the coemcients of Eq. (5) in which the condition £V ílj = 1 
is taken into account, and x is a set of N B — 1 independent, uniformly distributed random 
variables. The distribution of the quantities 5fi^ in this case is obtained from the expression 
P p = (í (/(«» -/f)>: 

pf = [ dkeW^) ri\ A ^ ( 7 ) 
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where every factor A p (x) = e p cos \ mp 1 ^ i determines the contribution to the integral 
from the corresponding periòdic orbit p. Thus, Eq. (7) gives the exact expression for the 
distribution Pj ^ in terms of the periòdic orbit theory. It is important to point out that 
the properties of the asymptotic distributions of trigonometric sums of form (6) are one of 
the traditional areas of research of mathematical statistics (see, e.g., [9,10] and references 
therein). In particular, it is known that separate terms (or groups of terms) of lacunary 
trigonometric series of form (6) can be considered as weakly dependent random variables, 
for which one can be establish a generalization of the central limit theorem, and consequently 
their sum is asymptotically Gauss distributed according to 

Pf = -U" 1 ^, (8) 

OV27T 

with the variance 



v 

The conclusion about the Gaussian form of the distribution of the fluctuations also appears 
to be applicable to this kind of spectral characteristics expansions of most of the regular 
quantum graphs (and other scaling systems), which are described by series of form (6) 
with constant coeficients. A hypothesis about the Gaussian nature of the distribution the 
spectral staircase fluctuations SN(k) = N(k) — N(k), confirmed by extensive numerical 
investigations, was previously proposed in [11] as the universal "central limit theorem for 
spectral fluctuations" applicable to general quantum chaotic systems. Owing to the existence 
of additional explicit expansions (5) this hypothesis, corroborated by the relation with the 
theory of weakly dependent random variables (trigonometric sums) can actually be extended 
to a much wider set of spectral characteristics. 

III. APPROXIMATE DESCRIPTION OF THE DISTRIBUTION FUNCTIONS 

Since the contributions of individual orbits to the series 5f^ behave as weakly depen- 
dent random variables, some physical simplifications are possible in Eq. (7). Expanding the 
exponentials A p (x), one can note that because expansion (5) is made in orthogonal harmon- 
ies, most integrals of the cross terms appearing from the product of the expansions A p (x) in 
Eq. (7). Contributions come only from the "resonant" terms for which one of the algebraic 
sums of the freqüències vanishes. The amplitude of these contributions decreases rapidly 
in the orders of the corresponding degrees of C p °\ that are proportional to the product of 
the corresponding number of scattering coefficients at graph vertices [6,7]. This argumenta- 
tion can be used to simplify the integral for P^- For example, in a simple approximation 
the contributions from resonances between different algebraically simple orbits can be dis- 
regarded. This is equivalent to untangling of the factors A p r(x) corresponding to different 
algebraically simple orbits, i.e., to the introducing an independent set of variables ay for 
every algebraically simple orbit. In this case, the distribution probability is represented in 
the form 
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Pj 0) = Jdke ik (^-^)\[Q pl (kC^)^ (10) 

p' 

where every factor 

Q p ,= f\ ik ^ ò ?'» co < vlfh v' x *' + ^)dx p ,, (11) 
Jo 

corresponds to the algebraically simple orbit p' and the sum with respect to v in Eq. (11) 
is calculated over orbits whose indices are múltiples of rh p . For a more crude description of 
the probability distribution profile, one can disregard the resonances between any distinct 
orbits, which is equivalent to the introduction of an independent phase x p for every orbit. 
Under this assumption, the integral in Eq. (7) is separated into independent integrals and, 
as a result, we arrive at the simple expression 

Pj°> = í dke<f W -^) JI J (*;Cf ) , (12) 

where Jq(x) is the zeroth Bessel function. Distributions of form (12) appear in communi- 
cation theory, for example, when analyzing the intensity of interfering telecommunication 
channels, the theory of wave propagation in random media, and other fields where stochastic 
signal models are used [12,13]. It is also worth noting that, in the approximation of inde- 
pendent random contributions, the conditions of the Lindeberg-Feller theorem and central 
limit theorem are satisfied, which establish the normal distribution law for the sum of in- 
dependent random variables. For spectral expansions (6) these conditions on the variances 
a p = (pp°^ /2 of individual contributions are satisfied due to the exponential increase in 
the number of periòdic orbits and the uniform exponential decrease of the magnitude of the 
coefficients Cp°\ As a result, in the approximation of independent random contributions, 
distribution (12) has the same Gaussian form (8), with the variance a 2 = ^2 p Cp°^ 2 /2 < oo 
as that predicted in [9,10] and [11] for the case of weakly dependent variables. Such descrip- 
tion is applicable to the statistical properties of various spectral characteristics of the regular 
graphs beginning with their harmònic expansions [5-7]. For example one can consider the 
fluctuations ón^ = ^ [k n — fc n ) , of levels around the average value, which have form (5) 
with W> = 0, y4°' ) — — f > an d the coefficients 
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or the difference sfj„ = k n+m — k n with Sm,n = j^ m i V 9 ^ = í ^r~ an d the coefScients 

^í sta (?MT)- (14) 

Knowing the distributions of these quantities, one can describe more complex objects such 
as the correlation function of fluctuations (àn^Sn+m^, autocorrelation function R 2 {x), and 
the form factor K 2 (t), given by the expression 

^ = f E< e " íwr > = f E^^Sw, ( 15 ) 

-^o — -^o — 

m m 

where F^}(k) is the characteristic function of distributions of form (7), (10) or (12), which 
are obtained from expansion (5) for s mi „ with coemcients (14), and thus, 

oo 
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^) = r E*)' ( 16 ) 

m=l 

It is important that all above distributions are closed expressions consistently describing the 
spectral characteristics in terms of periòdic orbit theory. 



IV. SPECTRAL HIERARCHY 



As mentioned above, in general quantum graphs are not regular and so for them the 
spectral expansions of form (5) cannot be obtained directly. A generalization to the irreg- 
ular case can be obtained by using the relationship between the two neighboring separator 
systems and kn~^ and by applying Eq. (4) to f(k) = k at the (j — l)th level of the 
hierarchy: 




Here, N^\k) corresponds to the spectral staircase of the sequence Uf 1 . Bootstrapping of the 
sequences by UP (or N^^k) by N®(k), see Fig. 1) means that jV^' -1 ) (kH^ = n. 

Substituting expansion (1) for Aü" 1 ) (jfe^) into Eq. (17), and using ïéjP in the form 
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= f (n + «>) , (18) 
we obtain the oscillating part of fc„ in the form 

e i} = - e <# _1) sin K ü " i)n + . ( i9 ) 

Here, the zeroth term 

= ± («> - e) - ± (ow - (er) , (20) 

the amplitudes, 

^-^^-V^'-e + l), (21) 
and phases cpp ^ = lu 

ü-i) ^OV^x-l) /2for every level j are functions of the fluctu- 

ations 5n ^ and 5^2 1 at the preceding hierarchy level. 

Similar expansions are easily obtained for other spectral characteristics, for example, for 
„ü-i) _ í.ü-1) _ lü-i). 

= + f e ^ c - (« - f^- i} ) , (22) 
° ? 

with the zeroth term 

fü-i) _ jj) + f s ü) _ oü) A x 

./s — ò n,m ^ l ò n,ra b n,rn-l) A 

(vrm/Lo - ( s g> m + *£ 1>m )/2) - (^L - s£i, m ) , (23) 

where # = (5n } + <#2i)/2 and the expansion coefficients Dp,m^ are obtained from the 
corresponding expansion for Sn]m- The equations relating the neighboring sequences can 
also be considered as describing the transition of a single separating sequence fn^ from one 
hierarchy level to another. 

V. STATISTICAL DESCRIPTION OF SPECTRAL HIERARCHY 

As in the case of the regular graphs, the description of the stochastic properties of 
sequences such as 5rP or Sn]m is based on the observation that parsing through the indices 
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n in the arguments of harmònic functions (19) and (22) leads to the appearance of random 
variables x. The idea of finding the distribution functions for various spectral characteristics 
is based on using the structural relations between the separating sequences obtained above 
in order to relate the probability distributions at different hierarchy levels. Beginning 
with the distribution Py at the regular level, one can determine the distribution -P^ ^ at 
the next level and so on, ending with the last, physical level. 
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FIG. 2. Distribution of variances at the odd levels of the spectral hierarchy for the four-vertex 
quantum graph with r = 7. The sòlid lines are the Gaussian approximations of the numerically 
calculated histograms. 

As an example, let us consider the behavior of the sequences 8n . For simplicity, we 
treat the fluctuations 8$ and <^f2i as independent random variables 8i and 8 2 distributed 
according to Pg . Correspondingly, one can write for the density Pg* (8) 

P$~ x) = j 5 (5 - ói j -V) P^P^d8 1 d8 2 dx. (24) 

Using Eq. (19) and representing the delta functional in exponential form, we obtain 

Pt 1 \t) = fdke M (j[Ag-^{x t S li S 2 )dx\ , (25) 

Here, the factors Ap (x, 81,82) correspond to the terms of expansion (19), which are now 
explicit functions of fluctuations at preceding hierarchy levels, and (*) n u) denotes averaging 
over these fluctuations with the weight 



(8 U 8 2 , k) = e -*tf- l) C'i*) Pp) (íl) P U) (j y . 



(26) 
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The expression (25) generalizes regular expansions (7), (10) and (12) for the single-level 
hierarchy to the general expressions for j > 0, averaged over the disorder at the preceding 
levels. 
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FIG. 3. Development of the probability distributions for the distances between the nearest neigh- 

(i) "(i) ï(i) 

bors Sn = kn — t = 3. The maximum distance between the nearest neighbors at the j = 

level in this case is s max = 8.68, for the regular celi size tt/Sq = 2.28. 

We note that the argument ation concerning the Gaussian distribution form in Section 
2 [9,10] can be directly applied to the distribution of ói^ at the regular level. However, as 
shown in Fig. 2, the distribution of 8Ü at higher levels j > is also Gaussian- like. For 
other spectral characteristics, for example, Sn (see Fig. 3), the sequence of transitions of 
form (25) can lead to asymmetric (non-Gaussian) distributions. 




VI. DISCUSSION 



The method proposed in [5] for solving the spectral problem is based on establishing the 
structural relationships between the sequence of physical levels k n and the regular sequence 
kn specified as an explicit function kn^ = k^ r+1 \n). For quantum graphs, the regular 
sequence is given by (2) and relation to k n is established through the system of auxiliary se- 
quences ïé£\ bootstrapping k n with k^- r+1 \n). The spectral hierarchy thus obtained consists 
of the system of sequences VïP and transition equations (17) from k^ to kn ^ . 

This approach allows not only the description of the evolution of base sequences kn^ 
from low to high hierarchy levels, but also the complete probability description of spectral 
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characteristics in the framework of periòdic orbit theory including those that are not directly 
described by the Gutzwiller formula. In this case, it is possible to follow the development 
of the scales of spectral fructuations, distributing disorder over the intermediate hierarchy 
levels, gradually passing from less to more disordered sequences. While the base sequence 
is maximally ordered, the amplitude of fluctuations in each next sequence Uf 1 increases as 
the index j decreases, i.e. with the approach to the physical spectrum [5]. The minimum 
number of auxiliary sequences kn^ necessary for bootstrapping /ci r+1 ' ) with kn^ defines to the 
complexity of the spectral problem with respect to the given bootstrapping method. 

The above relation between the properties of the series of expansions (19) and the prop- 
erties of weakly dependent random variables [9,10] reveals the physical origins of the uni- 
versality of the distributions of different spectral characteristics following from the limiting 
properties of the sums of such quantities. The existence of a sufficient number of transitions 
between hierarchy levels of irregular systems and, correspondingly, of averaging processes 
over random phases and disordered sequences in Eq. (25) leads not only to the Gaussian 
shape of the distribution of probabilities Pj ^ [as, e.g., for 5N(k) and, correspondingly for 
5n \ see [11] and Fig. 2), but also to the appearance of more complex (e.g., Wignerian, see 
[3] and Fig. 3) distributions. 

It is also important that determining the fluctuation probabilities in form (25) makes it 
possible not only to follow the appearance of general, universal statistical relations, but also 
to describe in detail the specific features of distributions Pj j \ which present the individual 
properties of each particular system. 

Work supported in part by the Sloan-Swartz Foundation. 
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